AD-A234  729  cumentat,on  page 


Form  Approved 
OMB  No  0704-0188 


:cileCion  :»  ^luOir^  suqqeit 

Suit*  5 204.  A/tington,  ,  A  ll 


1.  AGENCY  USE  ONLY  (Lt'ivt  OUnk)  I  2.  REPORT  OATE 


I  4.  TITLE  AND  SUBTITLE 


cxjl  '~YvV\e.  oo  ^eoo?^ 


3.  REPORT  TYPE  AND  OATES  COVERED 

FINAL  I  Jul88  to  31  Dec  90 _ 

S.  FUNDING  NUMBERS 
reW?  AF0SR-88-0256 

61I02F  2304/A5 


|  6.  AUTHOR(S) 


Mark  Pinksy  and  Volker  Wihstut; 


AFOSR-TR- 


0274 


7.  PERFORMING  ORGANIZATION  NAME'S;  A.NO  AQDRESS(ES) 
Northwester  Univeristv 
633  Clark  Street 
Evanston,  II  60208 


8.  PERFORMING  ORGA-Vt-l'iCN 
REPORT  NUMBER 


19.  SPONSORING- MONITORING  AGENCV  NAME^j  AND  AOORESS{ES) 


10.  SPONSORING  MONiTOR.NG 
AGENCY  REPORT  NUMBER 


•  Ins  AFB  DC  20332-6443 


AFOSR-88-0256 


11.  SUPPLEMENTARY  NOTES 


Vk;  APRi  G  153  i  v;  M 


IJj.  distribution  availability  statement 


>12b.\  DISTRIBUTION  CODE  ~  iji 


Approved  for  pah  lie  rciu  ijl  , 
uiulribution  unlimited. 


In The  joYnt^work  with  Volker  Wihstutz  of  the  University  of  North 
Carolina,  we  have  investigated  the  Lyapunov  stability  of  systems 


defined  by  a  system  of  differential  equations  with  a  stochastic 
driving  term,  which  may  be  either  white  noise  or  real  noise  .  In  the 


first  case  we  showed  that  for  nilpotent  systems  it  is  possible  to 
compute  an  arbitrary  number  of  terms  in  the  asymptotic  expansion  of 
Lyapunov  exponent  in  fractional  powers  of  the  noise  coefficient  when 
this  tends  to  zero.  This  includes  the  important  case  of  critically 
damped  oscillator,  which  had  not  been  treated  previously.  These 
results  were  then  extended  to  the  case  of  the  same  nilpotent  system 
driven  by  a  finite-state  Markov  noise  process.  This  was  obtained  by  a 
method  of  homogenization,  using  techniques  previously  established  to 
study  the  central  limit  theorem  for  functions  of  a  centered  Markov 
chain.  It  was  shown,  as  in  the  case  of  white  noise,  that  the  Lyapunov 
exponent  admits  an  expansion  in  fractional  powers  of  the  noise 
parameter,  and  that  the  first  term  of  this  expansion  agrees  exactly 
with  the  result  obtained  in  the  white  noise  case. 


14.  SUBJECT  TERMS 


IS.  NUMBER  OF  P*CE> 


17.  SECURITY  CLASSIFICATION  18,  SECURO  Y  CLASSIFICATION  19,  SECURITY  CLASSIFICATION  Z0,  L'MlTATlON  C»  RAv, 

WSSklFIED  »!»»  Wt’Mm  m. 


nrwninmi  kcekm 


Lyapunov  exponents  and  rotation  numbers  of 
linear  systems  with  real  noise 

Mark  Pinsky* 

Northwestern  University 

VOIKER  WlHSTUTZ* 

University  of  North  Carolina  at  Charlotte 


■IIS  GRA4I 

wnc  tab 

Uu&sncvaieefi 


*7 - - - 

Dlstrifcu 1 1 on / _ _ 

Availability  CoSas 


5  y,  i  tVV' J 


11  Research  supported  by  Air  Force  Office  of  Scientific  Research  TMst  >  j.ptciai 


1.  Introduction. 


A  large  literature  has  been  devoted  to  studying  the  asymptotic  properties  of  the  linear 
stochastic  system 

(1.1)  x;  =  AX,  +  eBX,F(G),  X0  =  x  €  Rd 

where  A,  B  are  constant  d  x  d  matrices  and  F(£,)  is  a  mean-zero  function  of  an  ergodic 
Markov  process  on  a  compact  state  space  M.  Of  particular  interest  is  the  top  Lyapunov 
exponent 

(1.2)  A(s)  =  limf"1  logSXti 

*fcc 


and  the  rotation  number,  suitably  defined  in  case  d  =  2.  To  analyze  this  system  it  is  noted 
that  the  joint  process  (X,.£t)  is  Markovian  on  the  product  space  Rd  x  M  with  infinitesimal 
generator 


(1.3)  L  =  L; v  =  G  +  AX  ■■  V  +  e(B X  -  V)F( 0 


Here  G  is  the  generator  of  the  noise  (0);  the  second  term  is  a  “systematic  derivative." 
making  no  reference  to  the  noise  process.  The  third  term  mixes  the  noise  and  state 
variables,  inviting  the  term  “noisy  derivative”  and  is  the  source  of  the  analytical  challenge. 
This  model  is  referred  to  as  a  “real-noise  driven  system.” 

It  is  well  known  that  the  noise  process  obeys  a  central  limit  theorem  in  the  form 


(1.4) 
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)ds  =$►  N{0.cr2t) 


(*10) 


a  normal  law7  with  mean  zero  and  variance  cr2t,  where  the  variance  parameter  a2  = 
— 2(6r-1i%  F)  and  the  inner  product  is  computed  in  terms  of  the  invariant  measure  of 
the  process  (£,).  Therefore  one  may  attempt  to  analyze  the  asymptotics  of  (A't,£<)  by 
studying  a  related  system  driven  by  white  noise 


(1.5)  d Yt  =  AYtdt  -f  sBY,  o  du-, 

where  (re<)  is  a  Wiener  process  with  mean  zero  and  variance  a2t.  The  diffusion  process 
(Vt)  is  Markovian  on  Rd  and  has  infinitesimal  generator 

(1.6)  I>  =  AY-V+  \s2(BY-X)2. 

One  may  conjecture  that  the  asymptotic  behavior  of  the  Lyapunov  exponent  for  the  real- 
noise  driven  process  (A"t)  is  equivalent  to  that  for  the  diffusion  process  [Yt),  at  least  to 
the  first  approximation.  It  is  our  purpose  to  carry  out  the  details  of  this  plan  in  cases  of 
interest. 
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2.  Nilpotent  Systems. 


In  a  previous  paper  [5]  we  investigated  the  Lyapunov  exponent  for  white  noise  systems 
with  a  nilpotent  deterministic  part:  Ad  =  0,  Ad~l  non- zero  and  B  generic.  This  includes 
the  free  particle  perturbed  by  multiplicative  white  noise  as  well  as  other  models  of  physical 
interest.  Recently  these  results  were  extended  to  systems  driven  by  “telegraphic  noise.’" 
where  M  =  {—1, 1},  by  Arnold  and  Kloeden  [7].  Now  we  can  show  that  these  results  can 
be  extended  to  the  general  real-noise  driven  system.  We  have  the  following  theorem. 

THEOREM  2.1.  Let  A(e),  r(e)  be  the  Lyapunov  exponent  and  rotation  number  of  the  2x2 
system  X't  =  [A  +  eBF(G)}Xt  where  A2  =  0,  A  non-zero  and  (£<)  is  a  finite-state  ergodic 
Markov  process.  Suppose  that  (Bex,e)  >  0  where  Ae  =  0,  e  non-zero,  (e,ex)  =  0.  When 
£  |  0  we  have 

(2.1)  A(£)  ~  Cj£2/3  r(£)  ~  C2£2/3 

for  positive  constants  C\,C2-  These  expansions  are  precisely  the  same  as  for  the  Lyapunov 
exponent  and  rotation  number  of  the  associated  diffusion  process  (Yt). 

We  compute  the  Lyapunov  exponent  by  the  “adjoint  method.”  This  consists  in  writing 
the  generator  in  terms  of  a  system  of  polar  coordinates  (p,  p>)  and  setting  Q(£,p>)  =  L(p). 
The  angular  process  (£<,«,?«)  is  ergodic  with  stationary  measure  N(df_  x  dp>)  and  the  usual 
formula  for  the  Lyapunov  exponent  is  \{e)  ~  fMxSd-i  Q(£,<p)N(d£  x  dip).  In  the  present 
case,  A(t)  may  be  characterized  as  the  unique  number  A  for  which  there  exists  a  function 
f{y,0  solving  the  equation  Lf  =  Q  —  A.  Indeed,  integrating  both  sides  against  N(d^  x 
d^)  shows  that  A  =  A(e).  Since  this  equation  may  be  difficult  to  solve  exactly,  we  may 
obtain  asymptotic  approximations  by  replacing  L  by  a  suitable  approximate  generator,  or 
equivalently  to  find  a  function  fe  and  a  number  Ae  such  that  Lfe  =  Q  —  Af  -f  0(Re)  for 
a  suitable  remainder  term  Rt.  Integrating  this  equation  against  A'(d£  x  dp>)  produces  the 
asymptotic  statement  Ae  =  A(s)  +  0(RC).  It  remains  to  find  the  approximations  fe.  A,. 

Re  •  _ 

To  do  this  we  apply  a  method  of  “homogenization.”  We  write  the  noisy  part  of  the 
generator  in  the  form  G  +  6V  and  show  that  this  is  approximated  by  a  diffusion  operator 
in  the  following  sense:  there  exists  a  second  order  operator  L0  :  C0C(R 2)  — ♦  C°°(R2)  and 
operators  L,  :  C°°{R 2)  — ♦  CX(R2  x  M)  ( i  =  1,2)  such  that  for  each  /  €  C°°{R2)  we  have 

(2.2)  (G  +  6V)(f  +  6L1f  +  62L2f)  =  62(LQf)  +  0(63 )  6  10) 

This  allows  us  to  reduce  to  the  case  of  a  white-noise  drive  system  for  which  we  know 
the  asymptotics  of  the  Lyapunov  exponent  and  rotation  number.  The  details  appear  in 
sections  4  and  5. 
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3.  The  Harmonic  Oscillator. 


We  now  compare  the  small-noise  behavior  of  the  stochastic  harmonic  oscillator  for 
the  cases  of  real  noise  and  white  noise.  The  real  noise  system  is  defined  by  X't  =  AXtdt  4- 

eBXtF(Zt)  where  ^=(_°7  o)andB=(l  o)  with  7  >  0.  We  assume  specifically 

that  (£t)  is  a  finite-state  Markov  process  with  self-adjoint  generator  (reversible  case)  and 
invariant  measure  i/(d£).  We  take  a  system  of  polar  coordinates  (p,  9)  with  11^/7  = 
ep  cos  y?,  X2  —  ep  sin  9.  The  stochastic  equations  take  the  form 

<f't  =  -n/7  +  £(^(6)/>/7]cosV« 

Pt  =  e[F Ut M  sin  ft  cos  <pt 

The  infinitesimal  generator  of  the  joint  motion  of  (£t,'rt>Pt)  is 

Lx  =  G  -  +  e[ F(£)/ /y]  ^cos2  +  sin  cos  y-— j  . 

The  white  noise  system  is  defined  by  the  Stratonovich  equation 
(3.2)  dYt  =  AYtdt  +  eBYt  0  dwt 


where  wt  is  a  Wiener  process  with  mean  zero  and  variance  cr2t.  Its  generator  is  given  by 
r  d  l  2  2,  2  &  • 

Iv  =  ~'/^  +  2£<’(cos  'i’a7  +  s,"'5cosv’a^) 

where  we  make  the  identification  a2  =  —2 (G~1F,F). 

To  obtain  the  asymptotic  form  of  the  Lyapunov  exponent  of  the  real-noise  driven 
system  we  look  for  an  approximate  solution  (/,  A)  of  the  equation  Lxf  =  Q  —  X  where 
QUi'p)  =  Lp  —  e[.F(£)/v/7]sinv?cos<)£’.  This  is  sought  in  the  form 

(3.3)  f  =  h+ef,+e'f, 

A  —  Aq  4  cAj  4  £  X2 


leading  to  the  conditions 

( G-  y/i-Q^)h  +  M?(0/v/T lc°s2  'rfo  =  )/ V^j  s*n -P cos Y  ~  Ai 

(G-SiA)h  +  im/^cos^f'  =  _x2. 
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This  leads  to  f0  =  0,  Aj  =  0,  /i  =  (G  -  -g^)  lF(£) an  v?cos^/v/7  which  may  be  solved  in 

ierms  of  eigenfunctions  ipk  where  Gxi'k  =  — for  k  =  1 . iV  with  xpi  =  1,  //]  =  0  and 

>  0  for  ifc  >  1  by  writing  r  (Osm^cos^  =  \  (F,  4>k)^k  sin  2^,  leading  to 


/.=(  iM)E 


t>i 


-  /zjt  sin  2-y  +  v/7cos2v? 


(F,vk)r{'k/(nl  +  47). 


Averaging  he  /,  equation  with  respect  to  the  normalized  measure  dtpv(d£)  and  integrating 
by  parts  leads  t  j 


-A2  =  (1/ x/7 )  J  F(Z)cos2  ipf[d<pv{d(,) 

=  -(2/v^)  J  F{^)sm'.pcos^f1d^v(di). 


From  the  above  spectral  representations,  we  can  compute  the  inner  product  of  fi  and 
F sin 2 y?  to  obtain 


/ 


fiF(Z)simpcos>fdfv{dZ)  =  ^  /^{F,  xpk)2(sm2  2^)/A{p\  +47) 

*>1 


with  the  result  ([4],  [6]) 

(3.4)  Ar/al  =  (l/47)E^<F,^)V(^  +  47). 

k>\ 

To  obtain  the  asymptotic  form  of  the  Lyapunov  exponent  of  the  -white-noise  driven 
system  v?  look  for  an  approximate  solution  (/.A)  of  the  equation  Lyf  =  Q  —  A  where 
Q(£,<p)  =  Lyp  =  (e2/27)cos2  ^cos2^.  This  is  sought  in  the  form 

f  —  fo  +  ?2h 
A  —  Aq  +  e2A2 


leading  to  the  conditions 


-V7/o  =  0 


1  d 

x/7/2  +  -a(cosv5^)2/o  =  ( ^2  /27)  cos2  <p  cos  2^  -  A2 
leading  to  the  choices  /0  =  0  and  the  well  known  result  (replacing  dy>  by  d^p/2z) 

A7hl,e  =  (<72/27)(2-)-1  [  cos2 v? cos 2^d^ 

(3.5)  ./-* 

=  (<t2/S7  )- 

If  we  make  the  identification  cr2  =  -2 (G-1F,  F),  then  this  may  be  written  as 

(3.6)  A^h“e  =  (l/47)^(F,v,)2/^. 


*>] 


These  computations  are  summarized  as 
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PROPOSITION  3.1.  We  always  have  the  inequality  A2eal  <  A2hlte. 

To  resolve  this  apparent  discrepancy,  it  suffices  to  consider  a  parametrized  family  of 
real-noise  processes  f  6~iF(^sf-s)ds.  When  S  J,  0  these  converge  to  a  Wiener  process 
with  mean  zero  and  variance  a2  =  —2(G~1F,F).  The  stochastic  equation  has  the  form 

Xl^AXt  +  elS-'FUit-WXt 

with  infinitesimal  generator 

Lx  =  AX  ■  V  +  6-2G  +  e6-lF(£)BX  •  V. 

To  obtain  the  corresponding  form  of  the  Lyapunov  exponent,  it  suffices  to  substitute  above, 
with  p.k  replaced  by  Hk&~2  and  F  replaced  by  F6-1,  Thus 

A?*'(«)  =  (i/41)r2^(/«r2)(F,fc)V[(^lrJ)2  +  47] 

k>  1 

=  (1/4-y)  \^nk(F,il'k)2/[nl  +  4764].: 

k>l 

When  6  J,  0  we  have 

Proposition  3.2,  limA2eal(6)  =  (1/4^)  ]T{F,V’k)2/W  =  A*1"10. 

610  k>  i 

Thus  we  retrieve  the  white  noise  result  in  the  CLT  limit. 


4.  Proof  of  Theorem  2.1  (special  case). 


We  are  given  an  ergodic  Markov  process  {£(f)}t>o  on  a  compact  state  space  M:  the 
infinitesimal  generator  is  denoted  G  and  the  invariant  measur®  v — thus  G*v  —  0  and 
G 1  =  0.  We  further  assume  that  the  Fredholm  alternative  is  satisfied  for  the  simple 
eigenvalue  zero,  i.e.,  the  inhomogeneous  equation  Gf  =  g  has  a  solution  provided  that 
Im  9(Ou(dt)  =  0;  the  solution  is  uniquely  determined  by  requiring  fMf(t Md£)  =  0. 
This  condition  is  satisfied  for  a  finite-state  Markov  process  or  for  Brownian  motion  on  a 
compacc  manifold,  for  example. 

Let  there  be  given  a  function  F(£)  with  mean  value  zero,  i.e.,  fM  F(£)i/{d£)  =  0.  Let 
(x(t).x'(t))  —  (xi(f),x2(f))  be  the  solution  of  the  second-order  system  x''(t)  =  £x(t)F(£(t)) 
with  the  initial  conditions  r(0)  =  xj,  x'(0)  =  x2.  This  is  a  Markov  process  on  the  product 
space  R 2  x  M  with  the  infinitesimal  generator 

0  A 

(4.1)  L  =  G  +  x2--  +  eF(OTiir-, 

OX  1  0X2 

The  top  Lyapunov  exponent  is  defined  by 

(4.2)  A  {?.)  =  limf1  log\/^i(02  +r2(/)2. 

(  J  OC 
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This  is  invariant  under  linear  change  of  coordinates  in  (xj,  x2)  space,  in  particular  the 
scaling  transformation  (xi,x2)  — ►  (xi ,  Cx2), 

We  introduce  a  system  of  “polar  coordinates’  by 


(4.3) 


xi=epcos<^,  ar2  =  Cep  sin  i^. 


We  make  the  identification  x(t)  =  xi(t),  x'(t)  =  x2(t)  and  consider  the  joint  process 
(Z{t),p(t),<p(t))t>o-  After  a  short  calculation  we  find  that 


f'(t)  =  -C  sin2  <p(t)  +  (e/C)  cos2  <p(t)F(((t)) 

p(t)  =  C  sin  f{t)  cos  f(t)  +  (e/C)  sin  <p(t)  cos  ip(t)F(£(t)). 


The  joint  process  (f(<),y?(<),  p(t))<> o  is  a  Markov  process  on  the  space  M  x  Rx  R  with 
the  infinitesimal  generator 


(4.5) 


L  =  G  +  eFU)xi~-  +  x2~ 

=  G  +  (— C sin2  +  (e/C)  cos2  ipF( f))J^ 

+  sin  v?  cos  <p(C  +  (e/C)F(£))-j£- 

ap 


•  in  fact  the  first  two  components  already  form  a  Markov  process,  but  we  shall  need  the 
full  generator  in  what  follows). 

We  write  the  genexotor  in  the  form  L  =  G  •¥  bV  +  D,  where 


^sin  ip 


d  2 

COSf-?—  +  COS  If 

dp 


A) 

d*J 


d  d 

D  =  C  sin  f  cos  f  - C  sin2  f-^- 

op  Of 


and  6  =  e/C.  We  refer  to  G  as  the  noise  generator ,  V  as  the  noisy  derivatives  and  D  as 
the  systematic  derivatives.  We  also  note,  for  further  reference,  the  function  Q{f,t)  defined 
as  Lp(  f,£)  is  computed  as 


(4.6)  Q{f,  £ )  =  C  sin  f  +  (e/C)F(f )  sin  f  cos  f. 

In  previous  approaches  to  stochastic  Lyapunov  stability,  one  solves  approximately  the 
equation  Lf  =  Q  -  A  for  suitable  /  =  f(f,0  &nd  A  €  R.  In  the  present  case  this  is 
not  directly  possible,  because  of  the  presence  of  the  noise  variables  F(£)  which  intervene 
both  in  the  generator  and  in  Q{f-0-  Therefore  we  apply  a  process  of  “homogenization"1 
to  replace  the  operator  G  +  6V  by  a  suitable  diffusion  operator  in  the  (f,p)  space  and 
ultimately  replace  Qif.O  by  a  function  Q{p)  related  to  a  suitable  diffusion  process. 
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PROPOSITION  4.1.  There  exists  a  second-order  differentia]  operator  Lq  in  the  (p.  p)  vari- 
a bles  with  the  following  property:  for  any  f  €  C°°(R  x  R)  there  exist  correctors  /,  € 
C°°(R  x  Rx  M)  (i  =  1,2)  such  that 


(4.7) 


[G  +  6V](f  +  6fj+S2f 2) 

a  ^  - 

a-  +  sin 
dp 

=  62(L0f)(p,p)  +  0(63)  (6  X  0) 


( 


G+  UF(0(cosV 


ip  COS  ip 


dp 


)] 


(/  +  6f\  +  ^2/2) 


where  L0f(p,p)  =  |  a2  ^sin^cos^^  +  cos<^^^  /  and  o2  =  —2(G  ]F,F)  >  0.  The 
term  0(6 3)  is  estimated  in  terms  of  the  C3  norm  of  f. 

PROOF:  We  choose  the  correctors  in  order  to  cancel  the  term  of  order  6  and  to  render  the 
coefficient  of  62  independent  of  £.  This  requires  that  we  have  the  equations 


(4.8) 


Gfi  -  cos2  pF(Z)^.  -sm<pcosipF(Z)j^  =0 


df 


(4.9)  Gf2  —  cos2  <pF(£)^-  —  sinpcospF(^)—-  =  fimction  of  (t p,p). 

O'P  Op 


The  first  of  these  is  satisfied  by  taking 

/i  =  -Bit) 


2  df  df 

cos  p  tt-  4"  sm  p  cos  p  — —  , 
dp  dp  J 


where  H(£)  is  the  solution  of  GH  =  —  F,  normalized  so  that  JM  ff(£)i/(d£)  =  0.  With  this 
choice  of  fi  we  substitute  in  the  equation  (4.9)  for  f2  and  average  with  respect  to  v(d£). 
The  Gf2  term  drops  out  find  the  right  side  of  the  equation  is  found  to  be 


(4.10) 


1  2f  2  d 

2  °  (C0S  % 


+  sin  p  cosp— 
dp 


f  =:  L.f 


where  c*  =  2  /„  F({)H(£M<1£). 

Finally  /2  is  determined  by  solving  the  indicated  equation  (4.9)  subject  to  the  nor¬ 
malization  JM  h(Z,p, p)v(d£,)  —  0.  This  is  possible,  since  Gf2  has  been  arranged  to  be 
perpendicular  to  the  null  space  of  G*,  completing  the  proof. 

If  e/C  |  0  we  may  restrict  attention  to  the  “approximate  generator”  (ejC)2L0f  — 
C  sin2  p ^  -I-  C  s\npcosp-g-p.  In  order  for  the  terms  to  balance  we  are  led  to  the  equation 
(e/C)2  =  C  or  C  =  e2^3  as  we  had  in  the  diffusion  case  [5].  More  precisely,  we  consider  the 
white  noise  system  dx  —  x2dt ,  dx2  =  exj  o  dw  where  {u>(f)  :  t  >  0}  is  a  WTiener  process 
with  mean  zero  and  variance  parameter  o 2t.  For  this  system  the  infinitesimal  generator 

is  Le  =  +  e2y-  ( xi~§p~^  '■  if  we  take  the  polar  coordinate  system  xi  =  epcos  p. 

x2  =  Cepsinv?,  we  obtain  the  angular  equation 

(4.11)  dp  =  -C  sin  p  cos p  +  (at/C)2  cos2  p  o  dw 
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with  the  Q  function  ( Qe  =  Lep) 

(4.12)  Qe{<p)  =  —C  sin  ip  cos  <p  +  (ae/C)2  cos2  p  cos  2p. 

From  our  previous  work  [5],  we  know  that  if  we  choose  C  =  e2/3,  then  Le  =  e2b3L\, 
where  L\  is  hypoelliptic  with  invariant  measure  p  and  the  Lyapunov  exponent  is  Af  = 

Qi(<P)M(dtp)  >  °- 

In  order  to  find  the  Lyapunov  exponent  of  the  real-noise  driven  system,  it  suffices  to 
find  f(p,0  and  A~  such  that  £/(vs£)  =  Q{<p,0  —  A~  +  0(e),  To  do  this,  we  proceed  in 
three  steps. 

Step  1.  Let  pi{p,£),  P2(<P,0  he  the  correctors  such  that  ' 

(4.13)  ( G  +  bV)(p  +  6pi  +  62p2 )  =  62LoP  +  0(b3) 

In  our  previous  paper  [5]  we  showed  that  the  operator  —  sin2  p-§~  +  L0  on  [— tt,  7r]  has 
an  invariant  measure  p  and  satisfies  the  Fredholm  alternative  for  the  simple  eigenvalue 
zero  and  that  Ax  =:  f*nQ{'p)p(dp)  >  0,  where  Q(p)  =:  sini^cos^  +  L0p  =  sin',?  cos  + 
a2  cos2  <p  cos  2p. 

Step  2.  Let  h  =  h(ip)  be  the  solution  of  the  equation 

.  2  d  _ 

-  sm  <p- r-  +  L 

d P 

normalized  so  that  f*n  h(p)p(dp)  =  0. 

Step  3.  Let  h h2(p,^)  be  the  correctors  defined  above  for  the  function  h,  i.e.. 
(G  +  6V)(h  +  Shj  4-  b2h2)  =  b2L0h  +  0(63). 

Proposition  4. 2.  With  the  above  notations,  we  let 

(4.15)  /(?,0  =  h(<p)  +  S(hi  -  +  P{h2  -  P2)(V,0 

where  S  =  e/C,  C  =  e2/,3,:  Then  Lf{p,£)  =  Q(p,£)  -  e2^3X\  4-  0(e)  in  particular  the 
Lyapunov  exponent  A(e)  =  e2^3X\  4-  0(e),  e  l  0. 

PROOF:  Recall  that  L  =  G  +  SV  +  D  where  G  is  the  noise  generator,  V  contains  the  noisy 
derivatives  and  D  contains  the  systematic  derivatives.  From  step  1,  we  have 

( G  4-  SV)(p  4-  6 pi  +  62  p2 )  =  b2o2  cos2  p  cos  2<p  4-  O(e) 

L(p  4-  6 pi  4-  b2 p2 )  =  C  sin  <p  cos  p  4-  b2 o2  cos2  p  cos  2 p  4-  O(e) 

=  e2^)  4- 0(e) 

(G  4-  bv)(h  4-  bhi  +  b2h2)  =  b2L0h  4-  0(e) 

L(h  4-  Sh\  b2h2)  =  —  C sin2  ip  h\p)  4*  S2 L0h  4-  0(e) 

=  e2/3(Q(<^)-A1)  +  0(e). 
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Subtracting  these  two  gives 

L(h  —  p  +  6(hi  —  pi)  +  b2(h2  —  P2))  —  — £2^3Ai  +  0(e). 

Adding  Q(p,£)  —  Lp  to  both  sides  gives  the  desired  result.  To  complete  the  proof,  we 
may  argue  directly  in  terms  of  martingales:  writing  the  above  equation  as  Lf  =  — £2/3Aj  + 
0(e),  we  have  that  Mt  =:  f(p(t),tp(t),£(t))  -  J0  Lf(p(s),<p(s),((s)ds  is  a  martingale.. 
But  Mt  =  o(t)  when  t  T  00;  therefore  we  may  divide  by  t  and  take  the  limit  obtaining 
limtjooP(0A  =  —  Ai£2/3  +  0(e).  We  have  the  following  theorem. 

THEOREM  2.1(a),  The  top  Lyapunov  exponent  satisfies 

A(e)  =  lim  t~1p(t)  =  e2^3(sin^cost^>  +  a2  cos2  tp cos2^)„  +  0(e)  (e  j  0), 

<Too 


By  the  same  method  we  can  compute  the  rotation  number,  defined  as 
r(e)  =  limtfoo<”V(0-  To  do  this  we  verify  the  following: 

Proposition  4.3,  Let  k(p)  be  the  solution  of 

-  sin2  <pk'(p)  +  L0k  =  L0tp  -  sin2  <p  -  (Lop)^  +  (sin2  <f)„. 

Let  ki(p),  k2(p)  be  the  correctors  defined  above  for  f  =  p  +  k(p).  Then  L(p  +  k(p)  + 
el/3h(p)  +  e2/3k2(p))  =  e2/3(L0p  +  sin2  p)^  +  0(e), 

Noting  that  k(p)  and  the  correctors  pi,  p2,  fci,  k2  are  periodic  functions  we  have  the 
following  result: 

THEOREM  2. ltB).  The  rotation  number  is  computed  as 

r(e)  -  lim  p(t)/i  =  e2^3  (sin  tp  cos  tp  +  a2  sin2  p)n  +  0(e)  (e  i  0). 

<Ioo  ' 


5.  Proof  of  Theorem  2.1  (general  case). 

We  now  generalize  the  set-up  of  the  previous  section  to  the  stochastic  system 
(5.1)  *'(<)  =  [A +  e*X«t))*]*(f) 

where  x(t)  =  (x}  (t),x2(t))  and  A^0isa2x2  matrix  with  A2  =  0.  Without  loss  of 
generality  we  may  take  a  basis  in  which  A  =  ^  0  )  S°  ^a‘  system  has  the  form 

x\ (t)  =  x2(t)  +  eF(£(t))(buxi(t)  +  b\2x2(t)) 
x'2(t)  =  eF(£(t))(b2iTi(t)  +  b22x2(t)).. 
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Taking  a  system  of  polar  coordinates  x\  =  epcosy?,  x2  =  Cep  sin  y?  we  find  that 


p'  =  -Csin2  y?  -  £siny?(&n  cos y?  4-  612C  sin  y?)F(f) 

+  (e/C)cosp(b2i  cosy?  4-  b22C  sin  y?)F(£) 
p'  =  Csin  y?  cos  ip  +  ecosp(bu  cosy?  4-  &i2Csiny?)F(£) 
+  (e/C)sinp(b2i  cosy?  4-  b22C  sin  p(F(£). 

Making  the  choice  C  =  £2/3,  the  generator  has  the  form 

d_ 
dp 

+  G  +  £1/3V!  +eV2  +  e5/3V3 


L  =  -£2'3 


/  -  2  d 

tsln 


4-  sin  p  cos  p 


where 


2  d  „ 

V\  =  62i  cos  y?F(^)—  4-  621  F({) sin y? cosy? -s- 


dp 


dp 


d  d 

v2  =  (622  -  &n)siny?cosy?F(£)  —  4-  (in  cos2  y?  4-  622  sin2  y?)F(f  )  — 


V3  =  —  &12  sin  y?—  4-  612  sin y> cosy? 
op 


dp 


This  has  the  same  structure  as  in  the  case  of  white  noise.  If  the  coefficient  b21  is  non-zero, 
we  may  apply  the  homogenization  procedure  described  above  to  obtain  the  approximate 
generator  as 

e2/3b22lL0f  -  £2/3  | sin2  y?^  4-  siny?cosy>  j  4-  0(e). 

Applying  the  adjoint  method  we  again  obtain  the  result  in  the  form  as  stated  in  Theorem 

2.1. 
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